A Quantum Algorithm for the Hamiltonian NAND Tree 
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Abstract 

We give a quantum algorithm for the binary NAND tree problem in the Hamiltonian oracle 
model. The algorithm uses a continuous time quantum walk with a run time proportional to 
y/N. We also show a lower bound of y/N for the NAND tree problem in the Hamiltonian oracle 
model. 
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1 Introduction 



The NAND trees in this paper are perfectly bifurcating trees with N leaves at the top and depth 
n = \og 2 (N). Each leaf is assigned a value of or 1 and the value of any other node is the NAND 
of the two connected nodes just above. The goal is to evaluate the value at the root of the tree. 
An example is shown in figure ([TJ) . Classically there is a randomized algorithm that succeeds after 
evaluating only (with high probability) iV 753 of the leaves. This algorithm is known to be best 
possible. See [I] and references there. 

leaves 




root 

Figure 1: A classical NAND tree. 



As far as we know, no quantum algorithm has been devised which improves on the classical 
query complexity. However there is a quantum lower bound of y/N calls to a quantum oracle [4]. 
In this paper we are not working in the usual quantum query model but rather with a Hamiltonian 
oracle [U [3] which encodes the NAND tree instance. We will present a quantum algorithm which 
evaluates the NAND tree in a running time proportional to y/N. We also prove a lower bound of 
y/N on the running time for any quantum algorithm in the Hamiltonian oracle model. 

Our quantum algorithm uses a continuous time quantum walk on a graph [5]. We start with a 
perfectly bifurcating tree of depth n and one additional node for each of the N leaves. To specify 
the input we connect some of these N pairs of nodes. A connection corresponds to an input value 
of 1 on a leaf in the classical NAND tree and the absence of a connection corresponds to a 0. See 
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the top of figure @. Next we attach a long line of nodes to the root of the tree. We call this 
long line the "runway" . See the bottom of figure ([2]) . The Hamiltonian for the continuous time 
quantum walk we use here is minus the adjacency matrix of the graph. As usual with continuous 
time quantum walks, nodes in the graph correspond to computational basis states. 




-MOM 



Figure 2: The full Hamiltonian Ho + Hjj 

We can decompose this Hamiltonian into an oracle, Ho, which is instance dependent and a 
driver, Ho, which is instance independent. Hd is minus the adjacency matrix of the perfectly 
bifurcating tree of depth n whose root is attached to node of the line of nodes running from — M 
to M. We will take M to be very large. See figure ([3]). 




leaves of bifurcating tree 
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Figure 3: The oracle independent driver Hamiltonian Hp 



Ho is minus the adjacency matrix of a graph consisting of the leaves of the bifurcating tree and 
the parallel set of N other nodes. Each leaf in the tree is connected or not to its corresponding 
node in the set above. See figure ([?]). The quantum problem is: Given the Hamiltonian oracle Ho, 
evaluate the NAND tree with the corresponding input. 

• ••••••• parallel nodes 

•••••••• leaves of bifurcating tree 



Figure 4: The Hamiltonian oracle Ho 



Our quantum algorithm evolves with the full Hamiltonian Ho + Hp, which is minus the adja- 
cency matrix of the full graph illustrated in figure ([2]). The initial state is a carefully chosen right 
moving packet of length L localized totally on the left side of the runway with the right edge of the 
packet at node 0. It will turn out that L is of order y/~N. We take M to be much larger than L, say 
of order L 2 . We now let the quantum system evolve and wait a time L/2 which is the time it would 
take this packet to move a distance L to the right if the tree were not present. We then measure 
the projector onto the subspace corresponding to the right side of the runway. If the quantum state 
is found on the right we evaluate the NAND tree to be 1 and if the quantum state is not on the 
right we evaluate the NAND tree to be 0. 
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We have chosen our right moving packet to be very narrowly peaked in energy around E = 0. 
(Note that E = is not the ground state but is the middle of the spectrum.) The narrowness of 
the packet in energy forces the packet to be long. If we did not attach the bifurcating tree at node 
0, the packet would just move to the right and we would hnd it on the right when we measure. The 
algorithm works because with the tree attached the transmission coefficient at E = is if the 
NAND tree evaluates to and the transmission coefficient at E = is 1 if the NAND tree evaluates 
to 1. The transmission coefficient is a rapidly changing function of E but for \E\ < l/(16yN) the 
transmission coefficient is not far from its value at E = 0. To guarantee that the packet consists 
mostly of energy eigenstates with their energies in this range, we take L to be of order y/N. This 
determines the y/N run time of the algorithm. 

Our algorithm uses the driver Hamiltonian Hv to evaluate the NAND tree. An arbitrary 
algorithm can add any instance independent Ha (t) to Ho and work in the associated Hilbert 
space. We will show that for any choice of instance independent Ha (t) the running time required 
to evaluate the NAND tree associated with Hq is of order at least y/N. 



2 Motion on the Runway 

Here we describe the evolution of a quantum state initially localized on the left side of the runway 
in figure ||2J| headed to the right. M is so large that we can take it to be infinite as is justified by 
the fact that the speed of propagation is bounded. First consider the infinite runway with integers 
r labeling the sites. The tree is attached at r = 0. We then have for all r not equal to 0: 

H \r) = -\r + 1) - \r- 1) r ^ 0. (2.1) 

For 9 > 0, e and e~ lrS correspond to the same energy (nonnormalized) eigenstate of equation 
(|2.ip with energy 

E(6) = -2 cos 9 (2.2) 

but the first is a right moving wave and the second is a left moving wave. We are interested in a 
packet, that is, a spatially finite superposition of energy eigenstates, which is incident from the left 
on the node and the attached tree. This packet will scatter back and also transmit to the right 
side of the runway. The packet is dominated by energy eigenstates, \E), of the form on the runway 

(r\E) = [e lr6 + R(E) e~ lre ] for r < 0, (2.3) 

(r\E) = T{E)e lr9 for r > 0. (2.4) 

(The states \E) do not vanish in the tree.) There are other energy eigenstates, but we will not need 
them. Furthermore from standard scattering theory we have 

(E(9)\E(9'))=2n8{9-9 r ). (2.5) 

Looking at r = we see that 

l + R(E) = T{E). (2.6) 
The transmission coefficient T(E) is determined by the structure of the tree. In particular let 

^> = <S|> <"> 

where "root" is the node immediately above r = on the runway, see figure @. Applying the 
Hamiltonian at \r = 0) gives 

H\r = 0) = -\r = -1) - \r = 1) - (root) (2.8) 
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and taking the inner product with \E) gives 

In the next section we will show how to calculate y(E) and show that if the NAND tree evaluates 
to 1 then y(0) = meaning that T(0) = 1 and if the NAND tree evaluates to then y(0) = oo and 
T(0) = 0. Unfortunately we cannot build a state with only E=0 since it would be infinitely long. 
Instead we build a finite packet that is long enough so that it is effectively a superposition of the 
states \E) with E close enough to that T(E) is close to T(0). We introduce two parameters e and 
D for which 

\T{E) -T(0)| < D\E\ for \E\ < e . (2.10) 

The parameters e and £> do depend on the size of the tree, but for the remainder of this section we 
only use e< 1 and sD <C 1. 

The initial state we consider, |^(0)), is given on the runway as 

(r|V(0)) = e tr7r/2 for - L + 1 < r < 

= for r < -L and r > (2.11) 
and vanishes in the tree. For this state 

(mmm) = o (2.i2) 

and 

(mw 2 \m) = - (2.i3) 

so the spread in energy about is of order 1/yL, However, for our purpose we will see that this 
state is effectively more narrowly peaked in energy, and (|2.13[) is really an artifact of its sharp edge. 
In fact most of the probability in energy is contained in a peak around of width X/L. Because we 
take 7r/2 and not —tt/2 in (|2.11[) we have a right moving packet with energies near 0. 
Evolving with the full Hamiltonian of the graph we have 

m))=e- iHt m0))' (2-14) 
We decompose |^>(i)) into two orthogonal parts, 

\m) = \Mt)) + \Mt)) (2.i5) 

where 

\Mt))= J ^e-^\E(9))(E(9)\m) (2.16) 

and \E(9)) is given in (|2.3p and (|2.4|) with the normalization (|2 . 5[) . The state \if>2(t)) is a super- 
position of the other eigenstates of H, that is, states incident from the left with < 9 < § — £ 
and ^ + e < 9 < it, states incident from the right and bound states with \E\ > 2 and (r\E) — * 
exponentially as r — * ±oo on the runway. We do not need the details of \ip2(t)) since we will show 
in the appendix that the norm of \ipi(t)) is close to 1 so the norm of \ip2 (t)) is near 0. In other 
words, \i^i(t)) is a very good approximation to \ip(t)). To ensure this we need 

L>-. (2.17) 

e 
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We would then expect that at late enough times we will see on the right a packet like the incident 
packet, but multiplied by T(0) and moving to the right, with the group velocity dE/d8 evaluated 
at 9 = tt/2 which is 2. In the appendix we show, if t is not too big, that for r > 

(»#(*)) = T(0) (r - 2^(0)) (2.18) 

plus small corrections. (In the appendix we also make sense of this equation for 2t not an integer.) 
To ensure (12.18[) we also need 

i>L> 2 e. (2.19) 

Since | -0(0)} is a normalized state localized between r = —L + 1 and r = 0, (|2.18[) implies that for 
t > L/2 

J2\(r\m)\ 2 = \T(0)\ 2 (2-20) 

plus small corrections so the cases T(0) = 1 and T(0) = can be distinguished by a measurement 
on the right at t > L/2. 



3 Evaluating the transmission coefficient near E=0. 



From the last section, equation (|2.9p . we see that we can find T(E) if we know y(E). We now show 
how the structure of the tree recursively determines y{E). Consider the tree in figure ^ ignoring 
the runway but including the runway node at r = 0. Except at the top or bottom of this tree every 
node is connected to two nodes above and one below. See figure (5) where a, b, c and d are the 
amplitudes of \E) at the corresponding nodes. Applying H at the middle node yields: 



from which we get 



E a = —b — c — d 
I 



Y 



E + Y 1 + Y" 

where Y = a/d, Y' = b/a and Y" = c/a. This is shown pictorially in figure (6). 

b • c 



(3-1) 
(3.2) 




Figure 5: The amplitudes at 4 nodes in the middle of the tree in an energy eigenstate. 




E+Y' +Y" 
Figure 6: The recursion for Y. 
The y we seek is Y at the bottom of the tree. To find it we recurse down from the top of the tree. At 
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the top of the tree there are 3 possibilities with the respective Y's obtained by applying the Hamil- 
tonian. The results are shown in figure (7). From figures (6) and (7) we see that Y(—E) = —Y(E) 
so we can restrict attention to E > 0. 



Y= 





Y= 



Figure 7: The values of Y's at the top of tree. 




We now show that the recursive formula for computing the Y's at E = is in fact the NAND 
rate. Looking at figure (7) and taking E — > positively we get figure (8). 



Y = 




Y = -oo 



Y = 




Y = 



Y = -oo 



Figure 8: The values of Y at the top of the tree as E 




Y = -oo 



We see that when we begin our recursion as E — > 0+ there are two initial values of Y which are 
-oo and 0. Returning to figure (6) with the different possibilities for Y'and Y" we get figure (9). 
Identifying Y = with the logical value 1 and Y = — oo with the logical value we see that figure 
(9) is a NAND gate. Accordingly the value of y(0) = Y(0) at the bottom of the tree (see 12 - T[) in 
figure ([2]) is the value of the NAND tree with the input specified at the top. 



Y = 




Y = -oo 



Y = 




Y = 



Y = -od 



Y = -oo 



Figure 9: The recursion for Y at E = implements the NAND gate. 




Y = -oo 



We now move away from E = to see how far from £ = 0we can go and still have the value of 
the NAND tree encoded in y(E) at the bottom. Near E = it is convenient to write Y(E) either 
as a(E) E or as —l/(b(E)E) and we will find bounds on how large a(E) and b(E) can become as 
we move down the tree. We observe that the coefficient of E in Y(E) (here we mean either a{E) 
or b(E)) on any edge is an increasing function of the coefficients on the two edges above, as long as 
\aE\, \a'E\, \bE\ and \b'E\ are not too big. To see this use the recursive formula in figure (6). There 
are three cases and we obtain figure (10) which makes the increasingness clear. 
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aE 



-l 

bE 



bE 

l-(l+o) bE 2 



aE 




-l 




1 



(1+a+a') £ 

Figure 10: The coefficient of E in Y is an increasing function of a, a', 6 and b' . 



-£ 2 



Now consider a piece of the bifurcating tree as depicted in figure (11). 




Figure 11: Y determined by four Vs above. 

There are 4 different Y's at the top of the piece, which determine Y at the bottom of the piece. 
Suppose for < E < 2~"/ 2 /16 each of the Y's at the top of the piece obeys 



0<Yi(E)<cE or 0< 



-1 
ME) 



< cE 



for some positive c with c + 1 < 4 • V where j < n/2. We will show 



< Y{E) < c* E or < 



Y(E) 



< c* E 



where 



(c* + 1) 



2(c+l) 



1 - 2""+ 2 V4' 

There are six cases to consider and here we exhibit the most dangerous case; see figure (12). 



(3.3) 

(3.4) 
(3.5) 




Figure 12: The most dangerous case. 

The reason we can use the coefficient c on the four top edges is the increasingness referred to earlier. 
Here Y(E) = d(E)E. Using the recursion from figure (10) we get 



d(E) = 



2c + 1 



< 



l-(l+ i- ( i + c e)eg (2c+l)£ 2 

2c + 1 

l-(l + 2c)(2c+l)£ 2 



(3.6) 
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using (1 + c) < 4 • 2"/ 2 and E 2 < ^ 2" 71 . So 

2c + 1 



d{E) < 



< 



1 - 4(1 + c) 2 E 2 

2(c+l) 
l-4(l + c) 2 £ 2 



2(c+l) 

^ l_2 2 i-"/4 ' 1 j 

using 1 + c < 4 • 2 J ' and E < 2~"/ 2 /16. From $£7) we get the left side of $£$$ with c* given in 
(|3.5[) . (From figure (10), the bounds on c and E also imply that the coefficients remain positive 
from one level to the next, so < Y(E).) The other cases yield smaller coefficients and we will not 
run through them here. 

At the top of the full tree Y{E) = E/(l - E 2 ) or Y(E) = -1/E as can be seen from figure (7). 
Given the restriction O < E < 2 _ ™ //2 /16 we can take 

ctop = 1 _ 1 2 _ n ■ (3.8) 
Moving down the tree we get after n/2 iterations of (|3.5p with j = 1, 2, . . . n/2, 

( C bottom + 1) = ( C top + 1) \TZ 2 2-n 7|J _ 2 4-n /4^ - " ' ^ 1 _ 2™-™ /ij ' ^' 9 '' 

Therefore 



2 n/2 

c k>ottom — ; _„/■? — ~ — (3.10) 

l-lE;4i2 2 ^-" 



(l + ^) 2™/ 2 



< A < 4 . 2 "/2 

1-1 

also justifying the assumption c < 4 • 2 J at each step. 

This means that \y(E)\ at the bottom is either less than |c bottom -E| or greater than | l/(c£ ottom .E) | 
for \E\ < 2~"/ 2 /16. The crucial VN arose because the coefficients of Y(E) could barely more than 
double as we moved two levels down the tree. 

Going back to the relation between y(E) and the transmission coefficient given in (|2.9p we 
summarize our results for this section: 



NAND = (reflect) \ y \ > \T\ < 8VN\E\ for \E\ < 
NAND = 1 (transmit) \y\ < As/N\E\ \T - 1| < 3y/N\E\ for \E\ < 



4 Putting it all together 

Here we combine the results of sections 2 and 3 and state the algorithm. First, the algorithm. Given 
the Hamiltonian oracle corresponding to an instance of the NAND tree problem we construct the 
full Hamiltonian Hp + Hjj which is minus the adjacency matrix of the graph in figure 2. We then 
build the initial state 

v r— — L+l 
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where r is on the runway. We choose 

L=jVN (4.2) 
with 7 an N independent constant and 7> 1. We let the state evolve for a time 

tiun = g ( 4 - 3 ) 
and then measure the projector P + onto the right side of the runway 



P + = J2\r)(r\- (4.4) 



If the measurement yields 1 we evaluate the NAND tree to be 1 and if the measurement yields 
we evaluate the NAND tree to be 0. 

According to the results stated in Section 2, the probability of getting a measurement result of 
1 is very near |T(0)| 2 if 

L>- and L > D 2 e. (4.5) 



From the table at the end of section 3 we can take 

1 



16ViV 

and 



(4.6) 



D = 8VN (4.7) 
and then the choice 7 > 1 ensures (|4.5I) . From the end of the appendix we see that the error 
probability of the algorithm is O (y^p (x) 1 ^) wmcn is O (\J~^j\ an d independent of N. 

By choosing 7 large enough we can make the success probability as close to 1 as desired. 



5 A lower bound for the Hamiltonian NAND tree problem 
via the Hamiltonian Parity problem 

Here we show that if the input to a NAND tree problem is given by the Hamiltonian oracle Ho 
described in section 1 then for an arbitrary driver Hamiltonian HA(t), evolution using Ho +Ha{£) 
cannot evaluate the NAND tree in a time of order less than y/~N. This means that our algorithm 
which takes time \/~N is optimal up to a constant. 

In the usual query model, the parity problem with y/~N variables can be embedded in a NAND 
tree with N leaves pQ. To see this first consider 2 variables, a and b, and the 4 leaf NAND tree 
given in figure (13) which evaluates to (1 + a + b) mod 2. Using this we see that with 4 variables 
a,b,c,d the NAND tree in figure (14) evaluates to (1 + a + b + c + d) mod 2. This clearly continues. 
Since we know that the parity problem for y/N variables cannot be solved with less than of order 
y/N quantum queries, we know that the NAND tree problem cannot be solved with fewer than of 
order y/N quantum queries. 




Figure 13: This NAND tree evaluates to (1 + a + b) mod 2. 
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Figure 14: This NAND Tree evaluates to (1 + a + b + c + d) mod 2. 



In our Hamiltonian oracle model the input in figure (14) becomes the oracle depicted in figure 
(15). Here the labels on the vertical edges on the top are always or —1. A label means that the 
edge is not included (just as in figure (4)) and accordingly the corresponding Hamiltonian matrix 
element is 0. A label —1 means that the edge is there (just as in figure (4)) and the corresponding 
matrix clement of the Hamiltonian oracle is — 1. 



a-1 



b-i 



c-l 



d-1 



b-1 



a-1 



d-1 



c-l 



Figure 15: The Hamiltonian Oracle for the NAND tree set up to evaluate (1 + a + b + c + d) mod 2. 
The coefficients are or —1 corresponding to the edge present or not. 

We see that the Hamiltonian parity problem can be embedded in the Hamiltonian NAND tree 
problem. We did this for the Hamiltonian NAND tree oracle considered in this paper, but it can 
also be done for the more general Hamiltonian NAND tree oracle described in the Conclusion. 

We will now prove a lower bound for the Hamiltonian parity problem, in a general setting (see 
also [3]), which can be used to obtain a y/~N lower bound for the Hamiltonian NAND tree problem. 
The oracle for the parity problem is a Hamiltonian of this form: 

K 

i=i 

The Hj operate on orthogonal subspaces Vj with 

Hj = Pj Hi Pj (5.2) 

where Pj is the projection onto Vj. Furthermore we assume ||-ffj|| < 1. For each j there are two 

possible operators Hj ai \ a,j = or 1. The string ax, ai ■ ■ ■ &k is the input to the parity problem to 
be solved. In the Hamiltonian oracle model, an algorithm can use (I5.ip . but has no other access to 
the string a. This is the most general form for the parity oracle Hamiltonian that we can imagine 
and it certainly includes the oracle we used to embed Hamiltonian parity in the Hamiltonian NAND 
tree. 

Choose an arbitrary driver Hamiltonian iJ^(i). Let S be an instance of parity, that is, a subset 
of 1, 2, ...K with a,j = 1 iff j is an element of S. Starting in an instance independent state |"0(O)) 
we evolve for time T according to the Schrodinger equation 

ij t \Mt))=Hs(t)\Mt)) (5.3) 

using the total Hamiltonian Hs(t) 

H s (t) = g(t) H OP + H A (t) (5.4) 

where \g(t)\ < 1. With the inclusion of the coefficient g(t) it is clear that the Hamiltonian parity 
oracle model includes the quantum query parity model. 
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Let |'0s(T)) be the state reached at time T. A successful algorithm for parity must have 



> 6 



(5.5) 



for some K independent S > if the parity of S and S' differ. We now show that this suffices to 
force T to be of order K. Our approach is an analogue to the analog analogue [T] of the BBBV 
method [2J. We write S ~ S' if S and S' differ by one element. Summing on (S, S") which differ 
by one element gives 



dt 



E II - |^M*)> = E 2Im vfe(*)l (H s (t)-H s , (t))\^(t)) 



= E 2Im(^(t)|±A J - (s ,s |^(t)) 



(5.6) 



where j (5, S") is the element by which S and S' differ, and Aj = g(t) (h^ - ijj 0) ) . The + sign 
means that j(S,S')eS whereas the — sign means that j(S, S') e S' . Now we have 



Jt E \\\Mt)) - liM*)> < 2 E \(Mt)\Pj(s, S ') &xs,s>) Pms')\^(t)) 



<4 E ||fys,sollte(t)>|| ' \\Pj(s,S')Ws'(t)) 



since ||Aj|| < 2. Using ab < (| a 2 + | fo 2 ) gives 



E ||i^(*))-i^(*))|| < 2 E 

For fixed S', j(S, S') runs over {1, 2, ... X} so 



\P j (s,s>)\il>s(t)) + P j{s ,sn\^s>(t)) 



E *W)iv&'(*)> 

s~s\S' fixed 



< 1 



and similarly for fixed S 1 and thus we have 



E 



< 2 • 2 



A" 



since there 2^ possibilities for S" (or £). 
We've shown that 

| E II l<M*)>-hM*)> 



and since 



we can integrate to obtain 



E ||rfe(o))-hMo)) 



E \\U's(T))-\^(T)) 



5~5' 



<4-2 



= 



< 4 • 2 A ' • T . 



A' 



(5.7) 



(5.8) 



(5.9) 



(5.10) 

(5.11) 
(5.12) 

(5.13) 



11 



For each S' there are K choices of S, so a successful algorithm requires 



E \\\MT))-Ws'(T)) 



and we have the desired bound 



T > KS/4. 



> 2 K ■ K ■ 5 



(5.14) 
(5.15) 



Conclusion 

We are not working in the quantum query model but rather in the quantum Hamiltonian oracle 
model. In this model the programmer is given a Hamiltonian oracle of the form 



N 



where the Hj's operate in orthogonal subspaces. Each Hj is one of two possible operators Hj 



with bj = or 1 and the string b\ ... bjy is the input to the classical NAND tree that is to be 
evaluated. The quantum programmer is allowed to evolve states using any Hamiltonian of the form 
g[t) Ho + HA(t) where the coefficient |<?(t)| < 1 and HA(t) is any instance independent Hamiltonian. 
The programmer has no other access to the string b. 

The algorithm presented in this paper uses a time independent Ho + Hp which is (minus) the 
adjacency matrix of a graph, so our algorithm is a continuous time quantum walk. We evaluate the 
NAND tree in time of order ^/N which is (up to a constant) the lower bound for this problem. 

After the first version of this paper appeared, a N 1 / 2+e algorithm in the quantum query model 
was found^ using our algorithm as a building block. As far as we can tell, our improvement leaves 
that result unchanged. 
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A Appendix 

Here we flesh out the claims made in section 2 and put bounds on the corrections to the stated 
results. Since we work with 6 close to ir/2 it is convenient to write 

9 = ip + n/2 (A.l) 

and 

E (tp + n/2) = 2siny> (A.2) 

so (|2.16[) becomes 

€ 

hM*)> = / ge- 2it8i ^|^+7r/2))(^ + 7r/2)|^(0)) (A.3) 

— £ 

where e is the small parameter introduced in (|2.10[) . Using (|2.3[) and (12. lip and letting r go to — r 
gives 



L-l 



(A.4) 



(£(^ + 7r/2)|V(0)> = -^ {e lrv + (-l) r R*(E)e-^}. 

V r— 

We now break \ipi(t)) into two pieces 

\M*)) = \M*)) + \Mt)) (A.5) 



where 



with 



and 



with 



I'M (/.)> = / ge-^ si ^A(^)|^ + 7r/2)) (A.6) 



L-l 



I ^— C - 1 e lLtp - 1 
= ~^Y e lr * = — — (A.7) 



i - 7> ! / ) = / J e" 2lt sin *> iT (J5) B(<p)\E(<p + tt/2)) (A., 



1 i-l 1 l-(-l) L e- lLv 
B(» = -= V (-l) r e-* rv = -= i ^- . (A.9) 

v r— 



Note that and IV'sW) are not orthogonal. Now 

1 



i%)| 2 < T \ x for M<£ (A.10) 
Lcos^ ±s 
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and \R(E)\ < 1 so 



is of order e/L and we have 



On the other hand 



Now 



W B {t)) 



o 



dtp 



(MWaM) = I g \a(<p)\ 2 . 



while 



2tt 



£P%)r + 1^)1 =- ^7 



1 /" 1 / sin ^Ltp 



2tt 



i \ sin \<p ) L 



< 



since sin f > £ and (siniiv 5 ) 2 < 1. Combining (|A.13|) . (|A.14[I and (|A.15|) gives 



1 ~ Te < (^(*)!^(*)) < 1 



so 



We now require 



L> - 



so the norm of I^aW) is close to 1 
Now 

llM*)) > I^a(*)) - llM*)) 



so 



= 1-0 



— \L 

Le'VL 



(A.ll) 

(A.12) 
(A.13) 

(A.14) 
(A.15) 

(A.16) 
(A.17) 
(A.18) 

(A.19) 
(A.20) 



where 0(a,b) means O of the larger of a and b. This justifies our claim in section 2 that \ipi{t)) is 
a very good approximation to the true evolving state \ip(t)}. Since 



\Mt)) + \Mt)) 



we also have 



\Mt)) 



o 



1/4 



(A.21) 



(A.22) 



We measure the state \ip(t)) on the right side of the runway, so we need a good estimate of 



r>0 



(A.23) 
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However we take t run to be near L/2, As long as Le 3 «1 we have 

oo 

]TM0| 2 = O(Le 5 ). (A.36) 

r=0 

To keep this small we need Le 5 <C 1, but this follows from our assumption that Le 3 <C 1 since 
£ <C 1. The assumption that Le 3 <§; 1 helps us to establish the translation property in (2.18) which 
simplifies the picture of what is going on. Also 

£ K(*)l 2 <\% l T (£) - ?X0)| 2 |v%>)| 2 . (A.37) 

r=0 _ r 

Using \T{E) - T(0)| < D\E\ for \E\ < e we get 

EKWI 2 < W^sin 2 ^^^ (A.38) 
^ J 2tt L sin (i<^) 



— £ 



L> 2 e 

= 0(— ) (A.39) 

from which we get our condition 

X > L> 2 £. (A.40) 

Now from (|2.11j) and (|A.7|) we have 

(r|^(0)) = i r / A{ip)e irv dip (A.41) 



So from (|A.27[) we see that 

a r (t) = T(0) (r - 2*1-0(0)) (A.42) 



but only when 2t is an integer. However if 2t = m + r with m an integer and < r < 1, 

,2 1 



— 7T 



- O (i) . (A.43) 

We have shown that for r > 0, (r|0(t)) is well approximated by (r\ipi (t)} which is well approximated 
by (r\tpA(t)} which is well approximated by a r (t) which has the form (IA.42I) so 12 . 18[) is justified. 
Furthermore, using (lA~25]) . (|A~28]) . (|A~33|) . (jA3U), (jA~39]l and fAiHl we have for t > L/2 

m)\P+\m) = \T(0)\ 2 + O (-±=, (^) 1/4 ) (A.44) 

which can be used to estimate the failure probability of the algorithm. 
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